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Abstract 

A Navier-Stokes computation is performed for a 
ducted-fan configuration with the goal of predicting 
rotor-stator noise generation without, having to re- 
sort to heuristic modeling. The calculated pressure 
field in the inlet region is decomposed into classical 
infinite-duct, modes, which are then used in either 
a hybrid finite-element /KirchhofF surface method or 
boundary integral equation method to calculate the 
far field noise. Comparisons with experimental data 
are presented, including rotor wake surveys and far 
field sound pressure levels for 2 blade passage fre- 
quency (BPF) tones. 

1 Introduction 

Tone noise in ducted fan engines, resulting from 
the interaction and response of the moving rotor- 
blade wakes with the stationary stator vanes, prop- 
agates both forward and aft. through the duct and 
radiates to the far field. This noise impacts both 
communities as well as passengers and crew flying 
in aircraft. Theoretical models for predicting rotor- 
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stator interaction modes have existed for quite some 
time (see, e.g.. Tyler and Sofrin 1 ), but most compu- 
tational techniques to date rely extensively on exper- 
imental measurements and analytical scaling tech- 
niques. Although such heuristic methods can pro- 
vide reasonable predictions of engine tone noise in 
many eases, they are far from perfect and it is diffi- 
cult to know where to turn when they are in error. 
Oil the other hand, predictions from first principles, 
although more expensive, hold the promise of offer- 
ing greater insight into and cont rol of the physical 
mechanisms behind the noise generation and propa- 
gation processes. As computers continue to become 
more powerful and inexpensive, the greater expense 
of the first-principles approach is becoming less of 
an impediment . 

Rumsey et a.L 2 first explored the use of the first- 
principles approach for predicting the one blade 
passage frequency (1 BPF) modes in an advanced 
ducted propeller (A DP) model. Only t he forward- 
propagating modes were considered. The Navier- 
Stokes equations were solved time-accuratelv on a 
moving grid with sliding patched interfaces between 
the rotor and stator grids. Turbulent flow was 
achieved through the use of a one-equation turbu- 
lence model with a wall function approach. The 
acoustic modes inside the engine were both gener- 
ated and propagated using this numerical method, 
and a particular mode of interest was shown to prop- 
agate well upstream of the rotor inside the engine 
with minimal attenuation. Separate acoustic codes 
were used to propagate these internal results out to 
the far field. 

The study of Rumsey et al. demonstrated the 
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feasibility of using a Navier-Stokes CFD code to 
bot h generate and propagate the tone noise within a 
ducted fan engine. However, the initial study was for 
l BPF modes, which are often stronger and hence 
easier to predict than higher harmonic modes. Also, 
most engines are designed with blade counts such 
that all 1 BPF modes are cut off (do not propagate), 
so in this sense the earlier ADP case was unrealistic. 

The current study applies the same methodology 
as Ref. 2 to a more realistic case for which the 1 BPF 
modes do not propagate but the 2 BPF modes do. 
Furthermore, more extensive experimental data is 
available for this new case, including rotor wake ve- 
locity profiles between the rotor and stator rows, al- 
lowing for a bet ter assessment of the first-principles 
method. This report details the progress and lessons 
learned to date, with the purpose of guiding future 
computational efforts in this area. Specifically, rotor 
wake profiles between the rotor and stator rows are 
compared with measured data, and forward- moving 
acoustic mode amplitudes within the duct are ana- 
lyzed. These in-duct acoustic results are then used 
to compute the far field noise in comparison with 
experimentally measured levels, using both a wave 
envelope/Kirchhoff method and a boundary integral 
equation method. 

2 Numerical Methods 

2.1 

Navier-Stokes method in interaction region 

The CFD code used in the current investigation 
is CFL.U), 3 a widely-used structured-grid upwind 
finite-volume method. It neglects viscous cross- 
derivative terms, which results in the thin-layer 
Navier-Stokes equations in specified coordinate di- 
rections. Third-order upwind- biased spatial differ- 
encing on the convective and pressure terms, and 
second-order differencing on the viscous terms are 
used; it is globally second-order spatially accurate. 
The CFL3D code can solve flow over multiple-zone 
grids that are connected in a one-to-one, patched, 
or overset manner, and can employ grid sequencing, 
multi-grid, and local time stepping when accelerat- 
ing convergence to steady state. Dpivind-biased spa- 
tial differencing is used for the inviscid terms, and 
flux limiting is used to obtain smooth solutions in 
the vicinity of shock waves, when present. No lim- 
iter was employed for this study. Viscous terms are 
centrally differenced. The flux difference-splitting 
(FDS) method of Roe is employed to obtain fluxes 
at l lie cell faces. 

The CFL3D code is advanced in time with 
an implicit three-factor approximate factorization 
method. The implicit derivatives are written as spa- 


tially first-order accurate, which results in block- 
tridiagonal inversions for each sweep. However, for 
solutions that use FDS the block-tridiagonal inver- 
sions are further simplified with a diagonal algo- 
rithm. Turbulence equations are solved uncoupled 
from the mean equations. 

The turbulence model used is the one-equation 
Spal art -Alim aras (SA) model,* 4 in combination with 
a wall function approach. The wall function ap- 
proach obtains a “pseudo"-wall viscosity that forces 
the law-of-tlie wall to hold. Details can be found 
in Ref. 3. The advantage of using the wall func- 
tion form is that the grid spacing does not need to 
be as tight near the solid walls (y+ ^ 8U as op- 
posed to t/+ % 1), so the stretching into the field is 
less severe. This is particularly important for acous- 
tic computations, for which good grid resolution is 
needed throughout the field to resolve the acoustic 
perturbations. 

A sliding patched interface connects the rotor 
and stator rows. At this interface, non-conservative 
interpolations are used, as described in Rumsey. 5 
In this reference, an engineering rule-of-thumb was 
developed for determining the maximum computa- 
tional time step allowable so that the acoustic waves 
that pass through the interface are not distorted. We 
adhere to this rule-of-thumb in the current work. 

For all the computations performed in this study, 
the CFL3D code was run time-accurately with three 
r-TS subiterations per time step performed in com- 
bination with multigrid, as described in Refs. 2 and 
3. 

2.2 Propagation technique to the far field 

In principle, the CFD code could be used to di- 
rectly calculate the far field noise. However, because 
of the code is only second-order accurate, such an ap- 
proach would require a prohibitive number of grid 
points. To make the problem tractable, a hybrid 
methodology is used that combines the modal data 
extracted from the CFD solution with more efficient 
means of computing the far field noise. 

Two approaches are considered. In the first, the 
modal data is input to Eversman’s Finite Element 
Model (FEM) duct radiation code. 1 ’ This velocity- 
potential based FEM code is then used to propagate 
the acoustic waves away from the inlet region; no at- 
tempt is made to include acoustic propagation from 
the exhaust flow. While the FEM code has been 
found to be accurate in the near field, phase infor- 
mation is not well represented in the far field. To 
overcome this difficulty, a Kirch hoff surface is con- 
structed using the near field FEM solution, and a 
Kirchhoff integral formulation is employed to com- 
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pute the noise at any point in the far field. This 
approach is described in detail in Ref. 7. This was 
also the approach used in the previous rotor-stator 
noise study." 

The second approach uses the ducted fan noise 
prediction code TBIEM3D. 8 9 Using boundary in- 
tegral equation techniques for solving the time har- 
monic convective wave equation, TBIEM3D calcu- 
lates the sound scattered by an infinitely thin, con- 
stant radius cylindrical duct in a uniform axial flow 
field. Fan noise generation is approximated by sim- 
ple spinning monopole or dipole point and/or line 
sources. In the present simulation, spinning line 
sources were used to generate the circumferential 
modes of interest . The line source strengths were 
adjusted to match as closely as possible the radial 
mode decomposition from the GFD. 

3 Results 

3.1 Description of the test case 

An advanced high bypass subsonic fan was de- 
signed and built by the Allison Engine Company 
under contract to NASA Glenn Research Center, as 
described by Woodward et al. U) The fan has 18 ro- 
tor blades with a diameter of 22 in. (0.559 in), and 
42 stator vanes. In the acoustic tests performed at 
NASA Glenn, various configurations of stator vanes 
were used; forward radial, aft radial, swept only, 
and swept plus leaned. For the present computa- 
tions, only the radial stators in the forward posi- 
tion are considered. (In the forward position, the 
stator leading edge is approximately 3.95 in. (0.10 
m) behind the rotor trailing edge.) The computa- 
tions use a fan speed of 5210 RPM, which is 50% 
of the nominal design speed (approach setting), and 
only forward-propagating acoustic modes are con- 
sidered. The current computations are run using 
a free stream Mach number of M — 0.2. This is 
higher than the M — 0.1 used in Ref. 10. However, 
the higher free stream Mach number primarily af- 
fects the flow external to the engine, including the 
stagnation point on the cowl. The flow inside the en- 
gine itself is primarily determined by the rotor RPM 
and is not expected to differ significantly for the two 
conditions. The Reynolds number is taken to be 
4.882 x 10° per m. 

With 18 rotors and 42 stators, Tyler-Sofrin infi- 
nite duct linear theory 1 can be used to determine 
the forward-moving modes that propagate. Duct 
acoustic modes are generally characterized as (m, n) 
modes by their circumferential and radial mode 
numbers m and n . As discussed in Ref. 1, the mode 
number in for a rotor-stator interaction can be ob- 
tained from m — hB -f ATT where h is the BPF 


harmonic, B is the number of rotor blades. F is the 
number of stat or vanes, and k is zero or any positive 
or negative integer. 

In the current study we are focusing on the 
upstream-moving 2 BPF in = —6 modes (h = 2 
and A* = -1). All 1 BPF modes are cut-off for this 
configuration. 11 The m = —6 modes are periodic 
over 1 / (3 1 h of the engine circumference (60°), and 
rotate in the direction opposite to that of the rotor. 
These are the lowest-order circumferential modes, 
and hence have the least stringent grid-spacing re- 
quirements for the GFD code: the higher the order 
the mode, the more points are required to accurately 
represent it. Given an average in-duct Mach number 
of M — 0.265 between the rotor and stator rows and 
M — 0.25 in front of the rotor, the only m ~ —6 
modes that propagate upstream without decay ac- 
cording to the infinite-duct theory are the first three 
radial modes. The axial wavelengths associated with 
each of these modes are: A(_r».i) = 0.087 — 0.089 m, 
A ( —o t i?) = 0.104-0. 106 m. and A ( _ G>3 ) = 0.151-0.157 
m. The first number in the range represents the 
wavelength associated with the higher Mach num- 
ber. Note that these wavelengths are theoretical val- 
ues for a constant Maeh-number infinite duct; there- 
fore they represent approximate levels expected in 
the engine, where the Mach number varies. 

3.2 Grid size and spacing 

For this case, it is necessary to have a grid that 
covers a full 60° of t he inside of t he engine. Thus, 3 
rotor passages and 7 stator passages must be mod- 
eled for the engine with 18 blades and 42 vanes. The 
grid for the current study is made up of these 10 
zones, plus 3 additional zones modeling the exterior 
of the engine, out to a far field approximately 37 
in. (0.95 m) in front, 62 in. (1.57 m) above, and 53 
in. (1.35 m) behind the engine cowl (the engine cowl 
is approximately 34.5 in. (0.88 m) long). An overall 
view of the far field grid is shown in Fig. 1. The 
grid size of each rotor passage is 321 x 57 x 41 (these 
numbers represent axial, radial, and circumferential 
spacing, respectively). The grid size of each stator 
passage is 129 x 57 x 33. There are approximately 
4.12 x 1(G total grid points in the 13 zones. 

To do a reasonable job in an engine acoustic com- 
putation, the GFD code has to adequately resolve 
two important aspects of the flow: the '"creation 
phase and the "‘propagation" phase. The "creation 
phase is the phase that occurs near the stator lead- 
ing edge, where the interaction of the rotor wakes 
with the vanes creates the duct acoustic modes in 
the first place. The "propagation'* phase is the ad- 
vection of these modes through the flow field. File 
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idea behind a first-principles approach is to use the 
Navier-Stokes equations to resolve both the creation 
and propagation phases within the highly nonlinear 
regime of the engine, where simpler linear acous- 
tic’ methods are not expected to be valid. Once the 
acoustic modes have been propagated away from the 
nonlinear regions, linear acoustic methods can then 
take over and be used to determine the noise in the 
far field. 

The grid size and spacing is a very crucial as- 
pect of an engine acoustic compulation. This point 
cannot be stressed enough. To apply the current 
first-principles approach to compute acoustic fea- 
tures, it is necessary to sufficient ly resolve the acous- 
tic* waves of interest as they propagate and the re- 
gion where the acoustic modes are initially gener- 
ated must be fine enough to adequately resolve their 
creation. Several close-ups of the grid are shown in 
Figs. 2 - 6. The grid was generated using TI(iER. 1J 
In the following paragraphs, we will point out some 
important elements of the grid and their relationship 
to the acoustic aspects of the computation. 

The grid is clustered near all solid walls, and 
stretches as it moves into the interior. As mentioned 
earlier, to avoid cells that are too large (because of 
t.oo-large stretching factors), we utilize a turbulence 
model wit h a wall function and use a minimum spac- 
ing near walls such that the level at the first grid 
point from the wall averages approximately 80. The 
current grid maintains fine axial spacing within the 
engine from the front of the inlet to past the leading 
edge on the stator. This axial spacing is fine enough 
so that there are at. least 25 - 30 points per axial 
wavelength for each acoustic mode of interest . The 
results for the model problem studied in Ref. 5 indi- 
cate that t his is an adequate resolution, at least for 1 
BPF modes. A view of the k — l planes of one rotor 
and one stator passage is given in Fig. 2. (Note that 
the figure distorts the inner cowl apparent shape, be- 
cause the grid lines sweep circumferentially to align 
with the rotor angled direction.) The axial spac- 
ing is necessarily clustered in the regions near the 
leading and trailing edges, but the average spacing 
elsewhere ranges from roughly 0.05 in. (0.0013 m) to 
0.12 in. (0.0031 m). Typical circumferential spacing 
in the rotor and stator zones between the two rows 
of blades is shown in Fig. 3. Clustering in the cir- 
cumferential direction is necessary both for the wake 
region (see below) as well as near the blades. 

In front of the engine inlet and behind the stators, 
t he Navier-Stokes code is no longer used to track the 
acoustic modes, so the grid is allowed to stretch ax- 
ially in these regions. This stretching accomplishes 
two objectives: it reduces the total grid size, and it 


diffuses potential spurious reflections from the exte- 
rior boundaries of the grid. (If one were to attempt 
this type of computation with only an unstretched 
internal-engine grid, the acoustic field would be con- 
taminated by reflections from t lie boundary unless 
special lion-reflective acoustic boundary conditions 
were employed.) 

Reasonably accurate resolution of the rotor wakes 
is necessary, because it is their interaction with the 
stator vanes that causes the tonal acoustic modes. 
To achieve high resolution, the rotor grids must be 
clustered near the approximate locat ion of the wakes 
as they convert aft. A view of the j — 35 planes 
of one rotor and one stator passage are shown in 
Fig. 4. The "wake cuts" of the rotor grid are ap- 
proximately aligned with the local rotor direction, 
so the grid lines follow the approximate paths of the 
wakes themselves. Additionally, the rotor grids ex- 
tend downstream to be wit hin 0.877 in. (0.022 m) of 
t he stator leading edge. As a result , the rotor wakes 
are in a “clustered-grid" region for a long distance 
before the wake passes through the rotor-stator in- 
terface, where clustering is minimized to improve 
inter-grid interpolations. 

Initial computations with what was originally con- 
sidered "reasonable" axial spacing near the stator 
leading edge led to results in which the (—6, 1) mode 
was not captured at all in the "creation" phase. It 
was subsequently determined that considerably finer 
spacing is required near the stator leading edge. A 
close-up of t he stator leading edge region in the cur- 
rent grid is shown in Fig. 5. It is still an open 
question whether even finer grid spacing in the re- 
gions surrounding the leading edge would signifi- 
cantly benefit the predict ion of the acoustic modes. 

In all of the computations done for this study, 
t he rotor tip gap is crudely modeled by collapsing 
(“zipping") the rotor grid zones for the upper 5 grid 
points. This zipping results in an average tip gap 
spacing of approximately 0.03 in. (7.15 x 1 0~ 4 m). 
which is similar to the actual tip gap measured in the 
experiment near the leading edge, but is over three 
t imes larger than the tip gap measured at the trail- 
ing edge. A view of the rotor tip is shown in Fig, (). 
looking aft . We originally attempted to model the 
actual tip gap spacings, but the finer grid spacing 
required near the trailing edge made grid generat ion 
more difficult and also caused excessive local stretch- 
ing factors that we decided were best to avoid. 

3.3 Computations 

Computations were performed at a nondinien- 
sional time step of At' — 0.05 (nondimensional- 
ized by characteristic length divided by free stream 
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speed of sound). The characteristic length in the 
CFD solution is 1 in. (0.0254 in) and the speed 
of sound is taken as c^ — 340.445 m/s. Dimen- 
sionally, this time step corresponds with a time 
step of A / = 3.73 x 10“'’ s. The frequency of all 
2 BPF modes is given by I\ = 2 BQ/c = 57.69 
in -1 , so the time corresponding to one period is 
1 p = 2i r/(c 1?c ,/v ) — 3.2 x 10 -4 s. (A given rotor blade 
rotates through 1 blade passage = 20° in 6.4 x 10 -4 
s.) Thus, the current time step yields approximately 
86 time steps per period for the 2 BPF modes. This 
number exceeds the 60 recommended in Ref. 5. 

As a first phase of the computation for this case, 
we investigate the ability of the CFD code to propa- 
gate a specified duct acoustic mode forward through 
the duct. We do this by removing the stator 
grid zones and specifying time- varying duct acoustic 
pressures at the exit face of the rotor grids, accord- 
ing to the formula: 

p( ?\ 9 < x) — AJ m ( A’ r r )exp/[c,^, I\ t — mO — I \ a •*’]• ( l ) 

A is the magnitude of the perturbation, J m is the 
Bessel function of the first kind and order m , K r is 
the radial wave number, and K a is the axial wave 
number. In essence, this computation is a check to 
test whether the grid resolution is sufficiently fine for 
the propagation part of the computation. In other 
words, given a duct acoustic mode of given strength, 
is the given grid fine enough to propagate it forward 
without significant attenuation to the rotor blades 
(where some of the mode is scattered into different 
frequencies 13 ); and, for the part of the mode that 
survives the passage through the rotor, is the given 
grid in front of the rotor fine enough to propagate it 
forward without significant attenuation to the duct 
inlet? 

As shown in Figs. 7 and 8, when an initial (—6, 1) 
mode acoustic amplitude of approximately .4 — 40 
Pa is imposed at the rear of the rotor grid, it propa- 
gates forward to the trailing edge of t he rotor blades 
with very little attenuation on either the fine or 
medium-level grid. Then, roughly a third of the 
wave strength makes it through the rotors. The 
wave is propagated forward on the fine grid all the 
way to the inlet with no attenuation, whereas on the 
medium-level grid it loses a small percentage of its 
strength. Therefore, the current fine grid appears 
to be fine enough for the propagation of the (—6, 1) 
mode of interest . 

Next, we run the full problem, including the sta- 
tors, and determine the ability of the CFD code to 
successfully generate the 2 BPF modes. As a first 
step, we compare the computed rotor wakes with 
LDY data taken at NASA Glenn Research Center. 


The first station corresponds to the LDY location 
closest to the rotor while the second station corre- 
sponds to the LDY’ location closest to the stator. 
The second station is located within the stator zones 
of the computational grid. Note that the LDY data 
was taken with the stator vanes in the aft position, 
whereas the CFD is performed with the stators in 
the forward position. However, the position of the 
stators is not expected to significantly impact the 
rotor wakes at the positions measured. Axial veloc- 
ity contours at constant axial positions are shown in 
Figs. 9 - 12, with computed contours corresponding 
to the fine grid. Generally speaking, the computed 
wakes are much narrower and exhibit more curva- 
ture than the LDY : measurements. Axial velocity 
defect profiles at a radius of approximately 8.5 in. 
(0.216 m) are shown in Figs. 13 and 14, correspond- 
ing to the dashed lines in Figs. 9 - 12. Although 
the computed width and general shape are different 
from the experiment. CFD does a reasonable job 
capturing the depth of the rotor wake as well as the 
overall trend of its evolut ion downstream. The axial 
velocity defects on the fine grid are not much differ- 
ent from the medium grid at the first station, but at 
the second station the wake depth of the fine grid is 
larger and in better agreement with experiment. 

We now' turn to an analysis of the duct acoustic 
modes generated by the Navier-Stokes code. This is 
done in a post-processing step by decomposing the 
instantaneous pressures output by CFL3D within 
the duct into its component duct mode strengths. 
The strength A of the (—6, 1 ) mode is shown between 
the rotor and stator rows in Fig. 15 and in front of 
the rotor row in Fig. 16. Axial distances in t hese fig- 
ures are given in terms of the distance forward of the 
rotor stacking axis. The mode strength approaches a 
quasi-constant level as it propagates forward toward 
the rotor trailing edge, and is roughly constant after 
it passes through the rotor row and travels forward 
through the duct. The wave-like structure is evident 
in both figures, and is close to the theoretical value 
of A fc* 0.09 m. 

Note in Fig. 15 that there is no distortion in the 
acoustic signal evident near the rotor-stator sliding- 
zone interface at x = —0.1256 m. This lack of dis- 
tortion indicates successful transfer of acoustic infor- 
mation across the patched interface. 5 Between the 
rotor and stator rows, the (—6, 1) strength from the 
fine grid averages roughly 20 Pa. and in front of the 
rotor row, it averages roughly 5 Pa. This 75% per- 
cent reduction in mode strength is similar to that 
seen in the earlier prescribed-mode case. There is a 
modest increase in magnitude of the mode strength 
between the medium and fine level grids in front 


5 

American Institute of Aeronautics and Astronautics 



of the rotor, indicating that further refinement be- 
yond the fine grid may still yield increased acoustic 
strength. Results for the (—(>,2) and (—6,3) modes 
in front of the rotor are shown in Figs, 17 and 18. 
The (—0.2) mode from the fine grid is very weak, 
with an average strength of only about 3 Pa. The 
(—0.3) mode has a higher amplitude of about 10 
Pa. The wavelengths of both the (—0,2) and (— 0,3) 
modes are close to the theoretical levels. Note that, 
among the three modes, the (—6, 3) mode shows the 
most significant, change in magnitude between the 
medium and fine level grids. 

Although not shown, as a further check on the 
acoustic computations, we studied the mode behav- 
ior at a given station in front of the rotor as a func- 
tion of time. (liven that t he mode behavior is kine- 
matically determined, we know theoretically how 
much it should rotate in a given number of itera- 
tions. Results confirm that the (—6, 1) and (—(>,3) 
modes rotate at approximately the correct rate and 
in the correct direction. The (—0,2) mode is too 
weak to yield meaningful analysis. 

4 Far Field Noise 

The final step in the current, study was to use the 
ext racted w — —0 modal strengths in the inlet duct 
forward of the fan (at j* = 0.1209 in) to predict the 
far field noise. The predicted far field directivity is 
shown in Fig. 19. Note that past approximately 00 
deg., the measured levels in the experiment are ele- 
vated due to aft-end noise, which was not included 
in the (TD computation. As can be seen, the com- 
puted sound pressure levels forward of 00 deg. are 
significantly lower than the experimentally measured 
levels, using both the F EM /Kirchhoff approach and 
t h e T B I E M 3 1 ) a pproacl i . ( Not e t hat TB IE M 3 D t ises 
an idealized duct, whereas FEM/KirchhofF uses the 
actual duct geometry. Thus, results are not expected 
to agree; for example. TBIEM3D yields higher SPLs 
at higher emission angles because of diffraction from 
the idealized duct’s sharp leading edge.) The cause 
of t he discrepancy between ( TD and experiment is 
therefore almost certainly related to the extracted 
modal data from the ( TT) solution - the computed 
source strengths in the inlet region are too low by 
about a factor of 3-6. 

As discussed previously, there is a 65%-75% reduc- 
tion in the modal amplitudes as the waves generated 
bv rotor-stator interaction travel forward through 
the rotor and into the inlet. Although not shown, 
if the transmission losses through the rotor are ig- 
nored, i.e.. if the modal strengths between the stator 
and the rotor (e.g., Fig. 15) are used instead of those 
in front of the rotor (e.g. Fig. 16), then the pre- 


dicted far field noise levels increase by nearly 10 dB 
and are much (‘loser to the experimental measure- 
ments. Other noise-prediction methodologies ignore 
the transmission losses entirely. 11 

Why the inlet modal strengths are underpredicted 
remains an open question. It may be that, the inter- 
action of the rotor wakes with the stators is not pre- 
dicted adequately because the predicted rotor wakes 
are much narrower than the measured wakes. This 
underpredict ion of wake width is possibly due to ei- 
ther inadequate physics in the turbulence model or 
under-resolution of the rotor wake vortical structure 
and its inherent unsteadiness. 

Thinking in terms of a Fourier transform of the 
rotor wakes in the circumferential direction, nar- 
rower wakes have a higher frequency content than 
wider wakes. Perhaps, then, the computed noise 
sources are biased towards harmonics higher than 
2 BPF. Alternatively, the losses that occur as the 
sound propagates forward through the rotors may 
be exacerbated by the numerical scheme. However, 
solutions obtained for two grid resolutions did not 
show a significant variation of the modal amplitudes 
with mesh size. 

5 Summary 

Performing acoustic computations on full rotor- 
stator ducted fan configurations using the Navier- 
St.okes equations is by no means a simple push- 
button task. The work presented here is the sec- 
ond of a series of computations performed to as- 
sess whether Navier-Stokes solvers can be used to 
directly capture the sound generation mechanisms 
of rotor-stator interaction, without resorting to more 
heuristic modeling techniques. The earlier study, for 
a configuration exhibiting 1 BPF tones, had a fair 
degree of success in predict ing the far field noise. In 
the current study, for which the lowest order propa- 
gating modes are 2 BPF, the predicted far field noise 
levels are about 10 12 dB low. It was determined 
that the propagation within the duct of the 2 BPF 
inodes was not as difficult as their cnation via the 
first- principles approach. The wake widths obt ained 
from the computation were considerably narrower 
than those measured in the experiment. This may 
be be the primary cause of the discrepancy in the 
computed and measured far field noise levels. Addi- 
tional research is needed to explore in greater depth 
the grid and turbulence model requirements in the 
interaction region where t he rotor wakes convert, and 
the acoustic modes are created. 
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Figure 1. Far field view of grid (k = 1 plane). 
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rotor-stator interface 
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Figure 2. View of rotor and stator grid zones (k = 1 
plane). 



Figure 4. View of rotor and stator grid zones (j — 35 
plane). 



Figure 3. View of rotor and stator grid zones ( i = 
250 plane in rotor zones and i — 10 plane in stator 
zones). 


Figure 5. Close-up view near stator leading edge 
(j - 35 plane). 
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Perturbation pressure (Pa) 




Figure 6. Close-up view of rotor tip region. 



Stator L.E. Axial distance (m) Rotor T.E. 



Fan Face Axial distance (m) Inlet 

Figure 8. Strength of forward-moving imposed 
(— (j, 1) mode in front of rotor row. 


Figure 7. Strength of forward-moving imposed 
(—6.1) mode between imposed location and rotor 
trailing edge. 
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Comnuted Axial Velocity Contours Ax» a i vetoc* 



Computed Axial Velocity Contours 
x = 4.99 5210RPM 



Figure 9. Computed axial velocity contours 1.2 in. 
(0.030 m) downstream of rotor tip trailing edge, 
within rotor zones; dashed lines show location for 
wake defect plots. 


Figure 11. Computed axial velocity contours 3.4 
in. (0.08(5 m) downstream of rotor tip trailing edge, 
within stator zones; dashed lines show location for 
wake defect plots. 


Measured Axial Velocity Contours Axial veioat 



Measured Axial Velocity Contours 
x = 4.99 5210 RPM 



Figure 10. Measured axial velocity contours 1.2 in. 
(0.030 m) downstream of rotor tip trailing edge, 
within rotor zones; dashed lines show location for 
wake defect plots. 


Figure 12. Measured axial velocity contours 3.4 in. 
(0.08(5 m) downstream of rotor tip trailing edge, 
within stator zones; dashed lines show location for 
wake defect plots. 
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Normalized Axial Velocity Normalized Axial Velocity 
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Figure 13. Normalized axial velocity defect 1.2 in. 
(0.030 m) downstream of rotor tip trailing edge, 
within rotor zones, at radius near 8.5 in. (0.21(5 m). 
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Figure 14. Normalized axial velocity defect 3.4 in. 
(0.080 m) downstream of rotor tip trailing edge, 
within stator zones, at radius near 8.5 in. (0.210 



Stator L.E. Axial distance (m) Rotor T.E. 

Figure 15. Strength of forward-moving computed 
(-0. 1) mode between the stator row and rotor row. 



Fan Face Axial distance (m) Inlet 

Figure 10. Strength of forward-moving computed 
(— 0. 1) mode in front of the rotor row. 
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Figure 17. Strength of forward-moving romputed 
(— (3.2) mode in front of the rotor row. 




Figure It). Far field sound pressure levels; narrow- 
hand directivities for 2 BPF tone along a 2.24 m (88 
in.) sideline, — 0.1. 


Figure 18. Strength of forward- moving romputed 
(— (j,«5) mode in front, of the rotor row. 
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